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This study focused on the capability of NASA Tetrahedral Unstructured Software System’s 
CFD code USM3D capability to predict the interaction between a shock and supersonic 
plume flow. Previous studies, published in 2004, 2009 and 2013, investigated USM3D’s 
supersonic plume flow results versus historical experimental data. This current study builds 
on that research by utilizing the best practices from the early papers for properly capturing 
the plume flow and then adding a wedge acting as a shock generator. This computational 
study is in conjunction with experimental tests conducted at the Glenn Research Center I’xl’ 
Supersonic Wind Tunnel. The comparison of the computational and experimental data 
shows good agreement for location and strength of the shocks although there are vertical 
shifts between the data sets that may be do to the measurement technique. 
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I. Introduction 


A s part of the NASA Fundamental Aerodynamics Program High-Speed Project, evaluation and improvement of 
computational fluid dynamics (CFD) sonic boom prediction capability has been a high priority. From the 
NASA Boom Workshop of 2008^ to the more recent AIAA Boom Workshop of 2014^, numerous conditions and 
CFD codes have been examined and compared. As this work progressed, researchers began looking into different 
components of supersonic flight. One of the remaining questions is how does the engine’s plume effect the boom 
signature? As a first step to answering that question, a study was conducted^ to ensure that the unstructured CFD 
code (USM3D) could properly capture a jet’s supersonic plume. Building on the lessons learned from that study, 
research was then conducted both computationally and experimentally to predict and model the interaction of a 
shock with an engine’s plume. 

The experimental component of this study was conducted at the NASA Glenn Research Center 1-foot by 1-foot 
supersonic wind tunnel (GRC 1x1 SWT) in Cleveland, Ohio. Since this tunnel’s capability to test small-scale 
models was developed to validate results obtained for supersonic nozzle concepts"^, it was chosen as the facility to 
obtain fundamental data on nozzle plume and shock interaction effects. The testing conducted used a simple jet 
flow and simulated tail shock generated by a wedge shaped shock generator. Wind tunnel test results include 
Schlieren images of the nozzle exhaust shock structure and near-field static pressure profiles as collected from a 
static pressure probe. 

The computational study was conducted using NASA Langley’s USM3D CFD code. The study was conducted 
by modeling the test section of the GRC 1x1 SWT with the jet and wedge installed. Additional cases included jet 
alone and wedge alone in order to attempt to isolate the effects of the jet supports and quantify the effect of the jet 
on the wedge’s shock signature, respectively. 


II. Numerical Techniques 

This CFD study used the NASA Tetrahedral Unstructured Software System (TetrUSS) for all the computations. 
This CFD suite, created and maintained by NASA Langley Research Center, includes a geometry setup tool 
GridTool, an unstructured grid generation program called VGRID, a grid post-processor named POSTGRID, and the 
flow solver USM3D.^ 

a) Grid Generation 

VGRID is an interactive, or batch, tetrahedral unstructured grid generation program. The grids produced by 
VGRID are suitable for computing either Euler or Navier- Stokes flow solutions. The grid spacing is related to the 
strength of user-defined sources placed in the domain. The methodology is based on the Advancing-Front method 
(AFM)^ and the Advancing-Layers method (ALM).^ Both techniques are based on marching processes in which 
tetrahedral cells are grown on an initial triangular boundary mesh and gradually form in the field around the 
geometry. Once the advancing front process is completed in VGRID, an additional post-processing step is required 
using POSTGRID to close any open pockets and to improve grid quality like removing the sliver cells. 

In 2008, a new version of VGRID was introduced.^ Software updates to VGRID included growth rates, surface 
sources, and volume sources. The outer boundary sources were no longer required for grid growth rates. The outer 
boundary sources were replaced with a user specified growth rate of the grid that dictates how the grids grow 
outwards from the surface and their maximum size. Additionally, new volume source types became available in 
VGRID. In the older version of VGRID, only point and line sources were available. Now, users can define sources 
that are based on a sphere, a cylinder, or a cone. 

In order to capture the plume and other flow characteristics properly, it is very important to source the nozzle 
configurations adequately and appropriately. For example, the cell spacing, distribution, and source divergence 
angle needs to be correctly defined to capture the shear layer. Adequate cell density is also needed to accurately 
capture the core length and the shock structure. A past study ^ showed that five cylindrical sources for the nozzle and 
the field provided adequate grid density for capturing the jet plume. A series of line sources are also placed on the 
nozzle lip to enable a smooth start of the shear layer. The cylindrical sources were placed to ensure that the nozzle. 
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plume core and expansion of the shear layer were properly captured. Additional information about the sourcing 
method can be found in reference 3. 

b) CFD Code Information 

The USM3D code^ is a cell- centered, finite-volume Navier-Stokes flow solver that uses multiple flux splitting 
schemes including Roe’s flux- difference splitting scheme^^ to compute inviscid flux quantities across the faces of 
the tetrahedral cells. Several options for turbulent closure are available: the one-equation Spalart-Allmaras (S-A) 
model^^ (with and without a wall function), and several two-equation models, including Menter’s Shear Stress 
Transport (SST) model. The minmod limiter, which is typically used for supersonic conditions, was used during 
this study to ensure numerical stability. For the flow solution cases conducted in this study, either the S-A or the k-8 
turbulence model with Sarkar’s compressibility and pressure dilation was used. This variant of k-8 model was 
shown to work well for modeling supersonic jets in references 3 and 13. 


III. Experimental Test Information 

The experimental testing component of this study was conducted at the NASA Glenn 1-foot by 1-foot 
Supersonic Wind Tunnel. The test section for this tunnel is 12-inches by 12.2-inches by 53.25-inches long. The 
tunnel can operate at 10 discrete speeds ranging from Mach 1.3 to Mach 6 by using contoured nozzle blocks 
upstream of the test section. For Mach 2, the tunnel can obtain Reynolds number per unit length between 1.3 and 13 
million. 

a) Jet Flow Simulator 

The jet flow model consisted of three major parts: the upper plenum, the strut, and the lower plenum (as shown 
in Figure 1). The upper plenum design was 1-inch outer diameter, and had an overall length from the nosecone to 
the nozzle exit of 8.21 inches. The support strut was 0.35 inch in thickness and contained nine passages to supply 
airflow. Air supply lines were 0.25-inch diameter and supplied 90 psia of pressure each to achieve 20 psia at the 
nozzle inlet for test conditions designed to simulate the environment at 50,000 feet altitude. Air was supplied to the 
bottom of the strut through a lower plenum by the centralized 450 psig pressurized air system at the NASA Glenn 
Research Center. 



Figure 1. Wind tunnel jet flow simulator cut-away. 
Dimensions are in inches. 
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b) Shock Generators 

The 5 -degree shock generator wedge was fabricated from 17-4 stainless steel to provide a leading edge shock 
that would intersect with the nozzle plume. The wedge was built with a 5-degree wedge angle (10-degree inclusive), 
6-inch chord, 8-inch span, and was mounted on a diamond shaped support strut 4 inches high from the wind tunnel 
wall. The wedge possessed sharp leading and trailing edges with a radius of 0.005 inches. A second, smaller wedge 
was built similar to the previous wedge except with a smaller chord of 2.5 inches instead of 6 inches, and a smaller 
span of 3.6 inches instead of 8 inches. 

c) Data Collection 

For this test, the data was obtained at Mach 1.96 and the experimental pressure signatures were measured at 
one-diameter away from the nozzle centerline, collected on the side opposite the model support strut. Data were 
collected at this location due to space restrictions in a small wind tunnel, and is consistent with previous 
experiments. The static pressure probe had a capability to travel axially up to 8.0 inches to capture pressure 
profiles, with an x=0 chosen as the nozzle exit plane. The actuator position and probe starting location were 
changed during the test, dependent upon where the shock interaction was located. The tunnel set-up can be seen in 
figure 2. 



Dimensions are in inches. 

Two static pressure probes were built, a 10-degree cone standard probe and a design, based on Pinckney. The 
10-degree cone standard probe design, in figure 3a, had difficulty capturing shock waves. This difficulty was 
demonstrated by a rounding of the peaks where shock waves occur, caused by the high cone angle and the long 
distance between the probe tip and the static pressure measurement ports where the measurement was not in the 
same flow environment as the probe tip. The Pinckney probe design (P=10; (0=1.5), figure 3b, had less sensitivity to 
angle of attack while achieving similar relaxation of static pressure at the location of the measurement ports. 
However, static pressure data collected with the Pinckney probe demonstrated a constant offset in Ap/poo of -0.08 
when tested in an empty tunnel, suggesting that the Pinckney probe has an offset in the Ap/poo measurement (Figure 
4). These probes were mounted on a 3/8 -inch thick airfoil strut, allowing 5.625 inches of probe length between the 
static pressure holes and the airfoil support. For all experimental results presented, the Pinckney probe was used. 

Further information about this test can be found in reference 4. 


4 


>ro 


P/Poo 





i 


Diam. * J 




I 

I 

L--—'"'” p/pcx, = 0.99 


a) Standard probe 


Tangent cone 


curve section 

Tangent conic \ ^ ^ static holes in 

curve section \ / tangent cone section 



^Conical Tangent cone ’^Tube 


half-angle, ro 
b) Pinckney probe 

Figure 3. Pressure probes used at GRC 1x1 SWT. 
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Figure 4. Experimental results for the standard and Pinckney probes in an empty tunnel at GRC 1x1 SWT. 


IV. Computational Test Cases 

The computational grids were fully viscous and resolved to the walls. The first cell height was sized to obtain a 
y+ of 0.2. Solutions were run on the grids until the residual of the L-2 norm reduced by at least 3 orders of 
magnitude. The runs were conducted at a freestream Mach number of 1.96, a Reynolds number of 271,526 based on 
the diameter of the nozzle (1 in.), nozzle pressure ratios (NPR) of 6, 8, 10, 12 and 14, and a total temperature of 529 

R. 


a) Jet alone 

The first case examines the jet modeled by itself within the test section of GRC 1x1 SWT (Figure 5). The grid 
consisted of 37.4 million grid cells. This work was done with the strategy that the effect of the jet support could be 
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removed from further computational results. Figure 6 shows a Schlieren image from the tunnel when running with 
the jet alone. As shown, there are numerous shocks caused by the tunnel, the nosecone and strut of the jet, and then 
reflections off the test section walls. Since computationally, only the test section was modeled, there should be less 
interference from other shock sources. The results from all of the CFD runs are shown in section V. 



Figure 5. Jet modeled in GRC 1x1 SWT test section. 
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Figure 6. Schlieren image of tunnel shock structure (NPR=8). 


b) Large Wedge Alone 

This case examines the 6-inch wedge modeled by itself within the test section of GRC 1x1 SWT (Figure 7). 
This case was examined to baseline the effect of the wedge alone in supersonic flow. The grid had 26.6 million cells. 



Figure 7. Large wedge modeled in GRC SWT test section. 


c) Jet & Wedge 

This case examines the 6-inch wedge mounted along with the jet within the test section of the GRC 1x1 SWT 
(Figure 8). This was originally the main piece of the research since this configuration would directly address the 
goal of the study, i.e., determine the effect of a shock wave on a jet plume. However, due to the other shocks 
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produced by the tunnel and various supports within the tunnel, this set-up was not tested. The computational grid 
had 47 million cells. 



Figure 8. Jet and large wedge modeled in GRC 1x1 SWT test section. 


d) Jet & Wedge as tested 

In order to experimentally test the 6-inch wedge mounted along with the jet, the wedge support was moved 
forward and reversed. Consequently, after the test, a grid was generated to match the tested configuration (Figure 9) 
with 33.9 million cells. 



Figure 9. Jet and large wedge modeled as tested in GRC 1x1 SWT test section. 


e) Jet & Small Wedge 

After the first entry into the GRC 1x1 SWT, results showed that even with moving the wedge forward, the 
expansion from the middle of the wedge and the shock off the back of the wedge were intersected by other shocks 
from the tunnel, thereby contaminating the results (shown in the circled area within Figure 10). It was decided to try 
a smaller wedge (2. 5 -inches) in a second entry to the tunnel. This new smaller wedge (Figure 1 1 and Figure 12) was 
modeled as tested and was able to experimentally capture both the expansion waves and parts of the aft shock. The 
computational grid contained 52.3 million cells. 



Figure 10. Schlieren image of jet and large wedge Figure 11. Schlieren image of jet and small wedge 
in GRC 1x1 SWT test section (NPR 8). in GRC 1x1 SWT test section (NPR 8). 
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Figure 12. Jet and small wedge modeled in GRC 1x1 SWT test section (NPR 8). 


f) Jet, Small Wedge & Full Tunnel 

The complete GRC 1x1 SWT was modeled (contraction + test section + diffuser) to accurately predict the 
boundary layer thickness and obtain the pressure signature. The goal was to compare the pressure signatures 
obtained by simulating the full tunnel versus the test section alone. Figures 13a and 13b show the tunnel top and side 
views respectively. The sourcing scheme for the tunnel is shown in figure 13c. Both line and cylindrical sources 
were used. The jet sourcing is similar to the procedure described in section Ila earlier. As can be seen, cylindrical 
sources have been placed above and below the nozzle, which extend to the tunnel roof and floor to capture the shock 
reflections. A cylindrical source was also extended from the nose of the model to the start of the test section to 
capture the bow shock eminating from the nose of the model. The source spacings were adjusted to ensure adequate 
grid point resolution on the tunnel walls. The volume grid size was about 36.5 million cells. 



n 


a) 

Top view. 

r 


1 

b) 

Side view. 



c) Sourcing in test section region. 

Figure 13. Jet, small wedge and sourcing modeled in full GRC 1x1 SWT. 
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V. Comparison of Computational and Experimental Results 

In this section comparison of wind tunnel data and CFD results are presented for both small and large wedge 
shock generators. Test cases modeled were: (a) jet alone, (b) large wedge alone and (c) jet plus large wedge both as 
originally proposed and (d) as tested, (e) jet plus small wedge in the test section alone and (f) in the full tunnel. For 
all the test cases, simulations were conducted for a Mach number of 1.96 and NPRs of 6, 8, 10, 12 and 14. 

As discussed in the experimental section, there is a possibility that the experimental data could be off by a value 
of 0.08. Therefore, for all of the plots containing experimental data, a blue symbol represents the experimental data 
with 0.08 added to it. 

Figure 14 shows a comparison between the wind tunnel data and CFD results for the jet only case. The data 
shown are at an NPR of 8 and CFD with the k-8 turbulence model. The CFD results compare well with the 
experimental results with an offset of 0.08. The CFD was able to capture the magnitude of the initial shock. 



Figure 14. Computational k-s and experimental results for jet alone at a 
NPR of 8 in GRC 1x1 SWT test section. 

Figure 15 shows CFD results for the jet only, the wedge only and the jet plus large wedge using k-E turbulence 
model for NPR of 6, 8, 10, 12 and 14. The blue line is the signature for the wedge only case. The plume shifted the 
signature of the wedge forward in the axial direction and reduced the wedge signature amplitude by approximately 
0.05, as seen in red circled area of figure 15a. 

Figure 16 shows comparisons between wind tunnel data and CFD results for the jet and large wedge in the GRC 
1x1 SWT for NPR of 6, 8, 10, 12 and 14. These CFD results were obtained using the k-E turbulence model. The 
CFD captured the location and magnitude of the shock at an x value of approximately 4.5 inches. The CFD results 
did not capture the peak at an x of approximately 6.5 inches that was caused by a shock reflecting from the upper 
wall of the wind tunnel. The CFD computational grid was fine enough to capture the nozzle and the wedge shocks 
but was coarse in the region between the plume flow and the tunnel wall (red oval in Figure 17). Since the desire of 
this study was to look at shock-plume interactions, priority was not given to capture the reflected shocks. If further 
detail is desired a finer grid is needed to be able to model shocks reflecting from wind tunnel walls. 

Figure 18 shows the comparison between computational results (S-A turbulence model) and wind tunnel data for 
the small wedge and jet in the GRC 1x1 SWT. The CFD results (red curve) follow the pattern of the wind-tunnel 
data (symbols) and match the magnitude of the shock and the overall nature of the expansion waves. However, the 
CFD shock slopes are found to be smaller than the experimental signature shock slope and slightly shifted to the 
right (i.e., downstream). 

Figure 19 shows the computational and experimental results for the jet and the small wedge with the fully 
modeled wind tunnel. These CFD results were obtained using the S-A turbulence model. Overall, there is good 
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agreement between the CFD results and the wind tunnel data. The CFD results follow the same pattern as the wind 
tunnel data with CFD signature peaks being consistently higher than experimental data. By modeling the full wind 
tunnel, any 3 -dimensional effects coming from upstream of the tunnel are captured and appear to provide a more 
accurate solution. Modeling the test section alone has the disadvantage that it will not see the full boundary layer 
growth and any 3 -dimensional effects from upstream of the test section. 

Figures 20 and 21 compare the computed density gradient and the experimental Schlieren images of the jet alone 
and the jet with the large wedge, respectively, at an NPR of 8 in the GRC 1x1 SWT. Figure 22 shows a comparison 
between the experimental Schlieren images and the computed density gradient for both the test section only and full 
tunnel configurations CFD computations of the small wedge at an NPR of 8 in GRC 1x1 SWT. Good qualitative 
agreement between Schlieren images of experimental and CFD results were obtained. The computed density 
gradient captured all the features present in the experimental image although the shock angles from the nacelle lip 
differ and a more pronounced shear layer definition. 
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c) NPR 10 




e) NPR 14 

Figure 15. Computational k-s results for jet alone, wedge alone and jet and wedge modeled in GRC 1x1 SWT 

test section. 
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Figure 16. Computational k-s and experimental results for jet and large wedge as tested in GRC 1x1 SWT test 

section. 
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Figure 17. Mach contours and the computational grid for the as tested large wedge and jet in the GRC 1x1 SWT. 
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a) NPR 6 



□ EXP 
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Figure 18. Computational S-A and experimental results for jet and small wedge in GRC 1x1 SWT test section. 
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c) NPR 10 




Figure 19. Computational S-A and experimental results for jet and small wedge with fully modeled GRC 1x1 SWT. 
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a) Experimental Schlieren results. b) Computational density gradient results. 

Figure 20. Computational and experimental results from jet alone at a NPR of 8 in GRC 1x1 SWT. 



a) Experimental Schlieren results. b) Computational density gradient results. 

Figure 21. Computational and experimental results from jet and large wedge at a NPR of 8 in GRC 1x1 SWT. 



a) Experimental Schlieren results. b) Comp, density gradient results. c) Comp, density gradient results. 

Test section only modeled. Full wind tunnel modeled. 

Figure 22. Computational and experimental results from jet and small wedge at a NPR of 8 in GRC 1x1 SWT. 
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VI. Concluding Remarks 

In this study, the NASA Tetrahedral Unstructured Software System flow solver USM3D was used to simulate 
the interaction between a supersonic nozzle jet flow with the shocks emanating from a shock generator. The CFD 
simulation modeled a wind tunnel test conducted at the GRC 1x1 SWT. A small and a large wedge were used to 
generate shocks in both the experiment and the CFD study. Simulations were conducted at the test section Mach 
number of 1.96 and the jet NPRs of 6, 8, 10, 12 and 14. The two turbulence models used in the study: the two- 
equation k-e turbulence model with compressibility correction of Sarkar and the S-A turbulence model. 
Computations confirmed that USM3D could predict the location and strength of the generated shocks. Overall good 
agreement between CFD results and experiment was obtained with the CFD signature peaks being slightly higher. 
However, due to questions about the experimental data, the actual pressure ratio value of the measured shock was 
not confirmed. A wind tunnel test is currently being conducted to verify the previous experimental results. 
Additionally, due to the contamination of shocks from other sources within the tunnel, i.e., the goal of the study, the 
computational study only attempted to match experimental results from x=2 to approximately x=6. 

The computational study shows that engine plume flow effects the shock signature by moving it slightly forward 
and dampening the peak pressure ratio values of the shock and it was also noticed that the wedge shock bends the jet 
plume flow upwards. Although preliminary experimental results (not shown) agree with these results, further 
investigation is planned. 

Good qualitative agreement between Schlieren images of experimental and CFD results were obtained. The 
computed density gradient captured all the features present in the experimental image and a more pronounced shear 
layer definition. 
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